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Experimental spectra of the CERN/SPS experiments NA44 and NA49 are fitted 
while using four different equations of state of nuclear matter within a relativistic 
hydrodynamic framework. For the freeze-out temperatures, Tf = 139 MeV and 
Tf = 116 MeV, respectively, the corresponding freeze-out hypersurfaces and Bose- 
Einstein correlation functions for identical pion pairs are discussed. It is concluded, 
that the Bose-Einstein interferometry measures the relationship between the tem¬ 
perature and the energy density in the equation of state of nuclear matter at the 
late hadronic stage of the fireball expansion. It is necessary, to use the detailed 
detector acceptances in the calculations for the Bose-Einstein correlations. 


1 Introduction 

The equation of state (EOS) of nuclear matter at very high energy densities is 
heretofore unknown. Many observables have been proposed as a signature for 
a possible new state of nuclear matter, i.e., the quark-gluon plasma (QGP), 
which could be formed within the very hot and dense zones of nuclear matter 
- the fireballs - in relativistic heavy-ion collision experiments. Among the 
proposed signatures for a QGP are Bose-Einstein correlations (BEG), which 
are correlations of identical hadron pairs. BEG functions are sensitive to the 
space-time dynamics of the fireball and therefore should give clues about the 
EOS, which governs the evolution of those fireballs. 

Among the many scientific contributions on BEG (for a recentmverview, 
c/., the book by WeinerEl) only few of them (c/., e.g., Rischke et al. □) address 
the particular role of the EOS of nuclear matter in the intensity interferome¬ 
try of identical hadrons, i.e., BEG. It is the purpose of this paper to further 
investigate the interplay between the EOS and BEG. In the following, we shall 
consider a framework of analysis which is based on relativistic hydrodynamics, 
because this approach allows for an explicit use of an EOS. 

In order to address the issues raised above, we shall use the simulati^on code 
HYLANDER-G 0 with different equations of state and Gooper-FryeQ freeze- 
out. Experimental single particle momentum distributwns from I58A GeV 
Pb-|-Pb collisions, measured by the NA44ffl and NAdQlD Gollaborations, are 
used as fit criteria for the hydrodynamical (one-fluid-type) calculations. 
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Figure 1: Speed of sound, Cq, ratio of pressure and energy density, P/e, and temperature, 
T, as functions of e, for the equations of state EOS-I (solid lines), EOS-II (dashed lines), 
EOS-III (dotted lines), and EOS Ib (dashed-dotted lines), respectively. The dots in plot (b) 
correspond for each EOS to the starting values of P/e with respect to the achieved initial 
maximum energy density e^ at transverse position = 0, (dots correspond to = 139 
MeV whereas the open circle corresponds to T/ = 116 MeV). The squares indicate the final 
values of P/e at breakup energy densities, €f. The dots A, P, C in plot (c) indicate the 
relationship between the temperature and the energy density at the late hadronic stage of 
the fireball expansion. 
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After a discussion of the space-time features of the particular fireballs we 
shall compare BEC of identical pion pair^to experimental correlation func¬ 
tions, which have been measured by NA44B. 


2 The Equations of State of Nuclear Matter 


In the hydrodynamic approach, virtually anj^ type of EOS can be considered 
when solving the relativistic Euler-equationsEl. The coupled system of partial 
differential equations necessary that describe the dynamics of a relativistic fluid 
(with given initial conditions) contains an equation of state which we write in 
the form 

P(e,n) = c^{e,n)e. (1) 

In Eq. (|^, the quantities P, e, and n are the pressure, the energy density, and 
the baryon density, respectively. The proportionality constant, c^, is in general 
a function of e and n. In the following, we shall assume that in particular the 
n dependence is negligible for the energy regime under consideration. From 
only the knowledge of^e speed of sound, cl{e), one can then calculate c^(e) 
by solving the integral^ ^ 

c2(e) = - / c2(e0de'. (2) 

£ Jo 


Assupiing an adiabatic expansion, the temperature, T(e), can be calculated 

fromy 


T{e) 


To exp 


'r clW)de' ■ 
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where Tq = P(eo) has to be sppcified for an arbitrary value of Cq. 

The first equation of stateM, EOS-I, which we use in the following exhibits 
a phase transition to a quark-gluon plasma at a critical temperature Tc = 200 
MeV with a critical energy density Cc = 3.2 GeV/fm^. The second equation 
of stateEJ, EOS-II, is also a lattice QCD-based EOS which has recently become 
very popular in the field of relativistic heavy-ion physics. This equation of state 
includes a phase transition to a quark-qluon plasma at Tc = 160 MeV with a 
critical energy density Cc « 1.5 GeVj jm?. The third equation of state^ EOS- 
III, has been extracted from the microscopic transport model RQMDEj under 
the assumption of complete thermalization, and does not include a transition 
to a QGP. We obtain a fourth equation of state, EOS-Ib, by changing the 
relationship between tc and Tc in EOS-I to Tc{ec = 1-35 GeV/fw?) = 200 
MeV. 

In Fig. 1 the four equations of state are plotted in many different repre¬ 
sentations. The particular parametrizations to the speed of sound which have 
been used here are published in Schlei et al. 113 . 
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Figure 2: Isothermes for the relativistic fluids governed by EOS-I, EOS-II and EOS-III at 
r± = 0, respectively. The outer lines correspond to a temperature, T = 140 MeV, and each 
successively smaller curve represents a reduction in temperature by AT = 20 MeV. The 
lines c = 1 represent the light cone. 
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Figure 3: Isothermes for the relativistic fluids governed by EOS-I, and EOS-Ib at r±_ = 0, 
respectively. For EOS-I the lines (beginning with the most outer curves) correspond to 
temperatures, T = 116 MeV, 120 MeV, 140 MeV, 160 MeV ... etc. For EOS-Ib, the 
outer lines correspond to a temperature, T = 140 MeV, and each successively smaller curve 
represents a reduction in temperature by AT = 20 MeV. The lines c = 1 represent the light 
cone. 
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3 Fits to Experimental Spectra and Freeze-Out 

In the following, we shall discuss five scenarios: we compare four calculations 
using EOS-I, EOS-II, EOS-III, and EOS-Ib for the fixed freeze-out tempera¬ 
ture, Tf = 139 MeV, and one calculation using EOS-I for the fixed freeze-out 
temperature, Tf = 116 MeV. 

In particular, it is assumed that due to an experimental uncertainty for 
the centrality of the collision, only 90% of the total available energy and the 
total baryon number have been observed. It is then possible to find initial 
distributions for the four equations of state, such that one can reproduce the 
single inclusive momentum spectra of 158 AGeV Pb-l-Pb collisions. Eor a 
freeze-out temperature, Tf = 139 MeV, the initial conditions and a large 
number of final single inclusive momentum djatributions for various hadron 
species have been published in SchleLet al. I3’E3 in comparison to the data 
measured by the NA44l3fl and NA49I3 Collaborations. The so far published 
results refer to calculations using EOS-I, EOS-II, and EOS-III. Results showing 
the fits of single inclusive particle momentum spectra using EOS-Ib wittur/ = 
139 MeV and EOS-I with Tf = II6 MeV will be published elsewhere EdI. 

It should be stressed, that all here discussed calculations result in single 
particle momentum distributions, which describe the corresponding data to 
the same extent very well. Although EOS-II was found in the cala^tions 
of hadronic transverse mass spectra to be too soft (c/. Schlei et al. we 
shall use it here also for the calculation of Bose-Einstein correlation functions. 
Before we discuss BEG, we have to discuss briefly the thermal evolution of the 
various fireballs. 

Eigs. 2 and 3 show isothermes for the relativistic Pb-|-Pb fluids governed 
by the different equations of state until freeze-out has been reached. The 
calculation using EOS-I with Tf = 139 MeV leads to a fireball of a much 
shorter liftime than the calculations using EOS-II and EOS-III with Tf = 139 
MeV. This behaviour is caused by much smaller freeze-out energy densities, 
e/, in the calculations using EOS-II and EOS-III compared to the calculation 
using EOS-I. We have £/ = 0.292 [0.126 (0.130)] GeV/fm^ when using EOS-I 
[EOS-II (EOS-III)]. A fluid that undergoes adiabatic expansion needs more 
time to reach the smaller freeze-out energy densities. 

Since EOS-II and EOS-III yield similar lifetimes of the fireball, in the 
following we attempt to increase the lifetime of the system which is governed 
by EOS-I. This can be acchieved by (a) using a smaller freeze-out temperature, 
Tf = 116 MeV, or (b) by hardening the EOS, i.e., using EOS-lb instead of 
EOS-I (without changing Tf = 139 MeV). For the latter two cases, we obtain 
e/ = 0.127 GeV/fm? (c/. Fig. 1 (c) and Fig. 3). 
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Figure 4: Projections of BEC functions for pairs emerging from the low p± (horizontal 

and vertical) acceptance setting of the NA44 detector The data points are data taken 
by the NA44 Collaboration The solid [dashed [dotted [dashed-dotted]]] lines correspond 
to the calculations using EOS-I [EOS-II [EOS-III [EOS-Ib]]] with Tf = 139 MeV, and the 
double dotted-dashed lines correspond to the calculation using EOS-I with Tf = 116 MeV. 
The values of the inverse width parameters, Ri {i = /,o, s), have been obtained from the 
Gaussian Bertsch-Pratt parametrization (see text). 
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4 Bose-Einstein Correlations 


The lifetime of theJOreball is, e.g., reflected in Bose-Einstein correlations of 
identical pion pairs E3. Within the Bertsch-Pratt variablesll3, a very sensitive 
quantity - regardless of whether one is interested in transverse expansion EH 
and BEC, or in the role of resonance decay ES in BEC - has always been the 
longitudinal projection of the two-pion correlation function, C 2 {qiong)- Fig- 4 
shows projections of BEC functions for 7r“'"7r''' pairs emerging frormthe low p± 
(horizontal and vertical) acceptance setting of the^A44 detector □. The data 
points are data taken by the NA44 CollaborationEi. 

It must be stressed, that in the calculations it was necessary to fully sim¬ 
ulate the experimental detector acceptance. I have checked, that if one does 
not follow the experimental prescription in generating the correlation functions 
one could obtain errors for the inverse widths and the overall strengths of the 
correlation functions in the order of 25%. 

Consistent with the expectation, the calculations using EOS-II, EOS-III, 
and EOS-Ib give similar lifetimes and therefore sufficiently large longitudinally 
expanded fireballs (see Figs. 2,3). In-Jiarticular, these calculations give - in 
addition to the excellent descriptionElE3llil of hadronic single inclusive momen¬ 
tum spectra (except for EOS-II) - a very nice reproduction of the pionic NA44 
BEC data. It should be stressed here, that a freeze-out temperature T/ = 139 
MeV was sufficient to achieve this agreement. 

Also consistent with expectation is that the calculation using EOS-I with 
Tf = 139 MeV gave a too small longitudinally expanded fireball (see above, 
and Fig. 2). On the contrary, it is on first sight surprising that a reduction of 
the freeze-out temperature to Tf = 116 MeV does not lead to large enough 
longitudinal extension of the fireball in the calculation which uses EOS-I. The 
reason for this result is the following: a freeze-out temperature reduction leads 
to a larger lifetime of the direct fireball, but the relative fraction of present 
resonance decay contributions is reduced (by about 30%); i.e., the resonance 
halo is reduced in size. Hence, the total fireball, which is a superposition of the 
direct (or thermal) fireball and the resonance halo remains moreless unchanged 
in size. 

In conclusion, by inspecting Fig. 1(c) we can see that only those equations 
of state which go through point C in the figure reproduce the experimental data 
on Bose-Einstein correlations good enough. Changing the freeze-out temper¬ 
ature shows hardly any effect. Therefore, the measurements of Bose-Einstein 
correlations tell us about which relationship between temperature and energy 
density is neccessary for a valid choice of an equation of state in the calcula¬ 
tions. Unfortunately, from the above made considerations it must be noted. 
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that two-particle BEC seen by itself cannot be used as a tool to determine a 
possible phase-transition to a QGP, because the BEC show no sensitivity to 
the structure of the EOS. 
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